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ABSTRACT 

The  mixed-mode  loading  of  a  rubbery  particulate  composite  is  studied  experimentally.  Linear  fracture 
mechanics  concepts  are  used  to  determine  the  initiation  of  growth,  the  initial  growth  direction,  and  the  subsequent 
growth  rate  for  a  range  of  mode  mixities.  The  fracture  toughness  locus  is  determined  to  be  elliptical  with  the  mode 
II  toughness  being  lower  than  its  mode  I  counterpart.  The  initial  growth  directions  correlate  with  maximum  strain 
energy  density  theories.  The  crack  growth  rates  can  be  modeled  effectively  using  an  equivalent  mode  I  crack, 
(keywords:  A:  particle  reinforcement,  B:  fracture,  mixed-mode  loading) 

INTRODUCTION 

A  crack  may  form  in  a  rubbery  particulate  composite  during  manufacture,  handling,  or  storage  of  a 
component.  Such  a  crack  will  be  in  the  mode  I  loading  orientation  during  this  transient  state,  but  will  later 
encounter  mixed-mode  loading  (also  called  combined  loading)  when  the  actual  service  loads  are  applied. 
Mixed-mode  loading  also  is  encountered  when  an  intrinsic  defect  causes  a  crack  to  f  n  at  an  angle  to  the  nominal 
™ode  I  loading  direction. 

An  early  study  by  Erdogan  and  Sih  used  maximum  tangential  stress  theory  to  predict  the  initial  crack 
growth  angle  (also  called  the  kink  angle)  and  the  fracture  toughness. ?  2  The  fracture  toughness  for  mixed-mode 
loading  can  be  described  by  two  parameters,  most  often  the  mode  I  and  mode  II  fracture  toughnesses,  Kjc  and 
Kjjc,  although  magnitude  and  phase  angle  are  also  used.  For  mixed-mode  loading  of  isotropic  materials,  an 
elliptical  fracture  toughness  locus  has  often  been  used3  Wu  analyzed  orthotropic  materials  under  combined  loading 
and  fit  the  fracture  toughnesses  with  a  similar  empirical  curve  fit4*5  Spencer  and  Bamby  generalized  the  results  to 
anisotropic  materials^5’6  Two-parameter  failure  loci  based  on  microstructural  assumptions  were  developed  by  Evans 
and  Hutchinson  and  by  Shih  and  Suresh.7,8  Numerous  experimental  works  have  confirmed  all  these  empirical 
curve  fits,  though  few  have  dealt  with  rubbery  particulate  composites,2-6,9 
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*■  '  For  many  materials,  the  fracture  locus  formed  by  plotting  critical  combinations  of  Kjc  and  Kjjc  will  be 

elliptical  due  to  changes  in  micromechanisms  that  cause  differences  in  the  toughness  as  the  mixity  changes. 
However,  the  relationship  between  the  fracture  toughness  locus  and  related  microstructural  phenomena  are  not  well 
understood,10  Wu  examined  micromechanisms  in  fibrous  composites,  and  Evans  and  Hutchinson  postulated  a 
surface  roughness  effect  to  explain  increased  toughness  with  mode  II  loading,4,7  Researchers  have  also  suggested 
effects  due  to  nonlinear  material  properties  such  as  plasticity,7,8,11  These  relationships  have  merit  for  some  classes 
of  materials  (e.g.,  metals),  but  may  not  apply  to  rubbery  particulate  composites. 

The  direction  in  which  the  mixed-mode  crack  will  begin  to  grow  has  also  been  studied.  Researchers  have 
proposed  various  criteria  to  determine  the  kink  angle,  or  initial  crack  growth  direction,5,2,5,10  However,  these 
theories  may  not  completely  explain  the  behavior  of  rubbery  particulate  composites,  which  employ  a  rubbery  matrix 
material  to  provide  structure  for  the  rigid  particles.  Because  of  the  matrix  material,  the  composite  can  sustain  large 
deformations.  Several  different  types  of  embedded  particles  make  up  a  large  volume  Vaction  (70-80%)  of  the 
composite. 

The  actual  mixed-mode  crack  behavior  in  these  composites  can  be  characterized  by  answering  three 
questions: 

1 .  Under  what  conditions  will  a  crack  start  to  grow? 

2.  If  the  crack  grows,  what  path  will  it  take? 

3.  At  what  rate  will  this  growth  occur? 

The  first  question  addresses  conditions  for  initiation  of  growth,  the  second  addresses  the  kink  angle,  and 
the  third  deals  with  the  growth  rate.  All  of  these  phenomena  depend  on  the  mode  mixity  and  are  important  to  the 
service  life  prediction.  As  an  example,  the  initiation  of  crack  growth  under  different  combined  loads  will  depend 
upon  two  components  of  fracture  toughness,  Kjc  and  Kjjq  The  fracture  toughness  locus  can  be  determined  by  a 
series  of  experiments.  Researchers  can  then  use  a  numerical  method  such  as  the  fin:  e  element  method,  along  with 

approximations  of  the  service  loads  and  geometry,  to  determine  the  complex  stress  intensity  factor,  K  =Kj+iKjj. 

This  stress  intensity  factor  can  be  plotted  on  the  failure  locus  to  determine  if  the  given  crack  will  grow. 

The  kink  angle,  or  initial  crack  growth  direction,  is  also  important,  because  it  determines  the  crack  path 
during  the  initial  stages  of  growth.  If  a  defect  or  other  feature  (for  example,  an  interface  or  a  free  surface)  exists 
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near  the 'crack,  growth  toward  it  and  the  resulting  interaction  could  result  in  accelerated  crack  growth  and 
component  failure. 

The  last  significant  consideration  that  affects  the  service  life  prediction  is  the  crack  growth  rate.  In  many 
instances,  the  growing  crack  will  grow  slowly  enough  to  be  noncritical  throughout  the  service  life.  However,  the 
need  to  prevent  catastrophic  fracture  requires  that  the  growth  rate  be  quantified  before  the  part  can  be  approved  for 
use.  Consequently,  more  accurate  appraisals  of  the  growth  rate  will  result  in  the  abi  lity7  to  use  components  that  were 
previously  questionable. 

The  current  experimental  work  addresses  these  issues  in  a  useful  and  simple  manner.  The  use  of  a  similar 
procedure  for  other  particulate  composite  materials  will  give  the  investigator  the  ability  to  predict  the  initiation  of 
growth,  the  initial  growth  direction,  and  the  subsequent  crack  growth  rate  for  mixed-mode  cracks.  This  is  achieved 
with  linear  fracture  mechanics  concepts.  Despite  the  viscoelastic  nature  of  the  rubbery  composite,  the  linear 
approach  provides  predictive  capability  (for  a  given  nominal  strain  rate)  and  can  be  used  without  a  thorough 
knowledge  of  the  time  dependent  behavior,  which  would  require  more  elaborate  experimentation  and  analysis. 
EXPERIMENTAL  PROCEDURE 

Figure  1  shows  a  typical  specimen  used  here.  The  cracks  are  at  an  angle  >3  to  the  mode  I  loading  direction. 
For  a  crack  in  an  infinite  plate  of  a  linear  elastic  material,  in-plane  traction  boundary  conditions  results  in  stress 
intensity  factors  given  by  12: 


Kj  =  crjna  cos2  ft 
K„  =  cryfna  sin/?cos/? . 


(l) 


Here  Kj  and  Kjj  are  the  mode  I  and  mode  II  stress  intensity  factors,  er  is  the  remote  nominal  stress,  and  a  is 
the  crack  length.  This  formula  applies  to  an  infinite  plate  with  traction  boundary  conditions.  In  such  a  case,  the 
crack  orientation,  ft,  is  identical  to  the  phase  angle  of  the  complex  stress  intensity  factor,  defined  by  ‘/'below: 


K=K,  +iKu  =Ke™ 


(2) 
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Here  the  complex  stress  intensity  factor  has  real  and  imaginary  parts  K[  and  Kjjbnt  can  also  be 
characterized  by  the  magnitude,  K ,  and  phase  angle  W.  The  actual  specimens  were  of  finite  width  and  had 
displacement  boundary  conditions,  but  finite  element  analysis  revealed  that  the  phase  angles  of  K  for  these  finite 
width  specimens  were  also  approximately  equal  to  the  crack  orientation  angle  We  tested  two  specimens  for  each 
of  five  different  crack  orientations:  0,  15,  30, 45,  and  60  degrees. 

A  displacement  controlled  testing  machine  loaded  the  specimens  at  a  fixed  rate  of  5.08  mm/min.  A 
special  fixture  applied  uniform  vertical  displacements  to  the  upper  and  lower  edges  of  the  specimens.  During 
testing,  the  specimens  were  videotaped.  We  synchronized  the  images  with  the  load-displacement  data  to  determine 
the  time,  displacement,  and  load  at  which  crack  growth  began.  Before  crack  growth,  significant  crack  tip  blunting 
occurs  due  to  the  large  deformations  that  the  material  can  sustain.  The  actual  growth  occurs  by  the  damage  of 
material  near  the  crack  tip.  Under  increasing  loads,  the  particles  near  the  crack  tip  separate  from  the  rubbery  matrix, 
resulting  in  the  nucleation  of  voids  which  then  coalesce  with  the  macrocrack. 

The  determination  of  the  boundary  conditions  at  initiation  is  easily  accomplished  using  the  videotape 
equipment  and  the  testing  machine.  These  boundary  conditions  (edge  displacements  in  this  case)  can  be  applied  to 
a  finite  element  model  that  uses  constitutive  properties  determined  from  tensile  tests  (see  Figure  2  for  an  example 
of  a  mesh).  The  finite  element  results  include  the  J  integral,  determined  here  using  the  domain  integral.  The  J 
integral  and  the  magnitude  of  K  are  simply  related  through 


Ef  =  E  {plane  stress),  E '  = 


E 

1- v2 


{plane  strain)  4 


(3) 


The  phase  angle,  ^  can  be  evaluated  using  the  crack  flank  displacements.  The  opening  and  sliding 
relative  displacements  of  the  two  faces  are  determined  in  the  limit  as  r  ->0  and  related  to  the  phase  angle  of  K.  A 
slightly  different  method  uses  the  stresses  Oyy  and  Oxy  along  the  positive  x  axis.  Analysis  using  these  methods 
revealed  that  for  the  plate  specimens,  the  crack  orientation  angle,  /?,  and  the  phase  angle,  *Fy  are  about  equal,  so  that 


the  phase  angle  *F was  taken  to  be  identical  to  the  crack  orientation  angle  /3  for  the  remainder  of  this  work. 

However,  for  determination  of  the  phase  angle  in  an  actual  structure,  the  use  of  crack  flank  displacements,  stresses 
along  the  *  axis,  or  some  other  method  is  necessary. 

The  use  of  eqn  (3)  for  determining  the  magnitude  of  K  and  the  knowledge  that  ft  =  IPprovide  a  way  to 

calculate  K  completely.  Using  the  boundary  conditions  (determined  from  the  videotape  and  load-displacement 
data)  at  the  point  of  crack  growth  initiation  and  incorporating  these  boundary  conditions  into  the  finite  element 
model  gives  the  fracture  toughness  in  terms  of  K\c  and  Kjjc,  or,  alternatively,  Kc  and  W 

Finite  element  models  were  used  with  experimental  results  to  determine  the  initiation  of  growth  in  the 
various  specimens,  but  the  kink  angles  were  determined  with  experimental  data  only.  Figure  3  shows  a  crack  in  a 
specimen  (the  crack  orientation  angle  was  60°).  This  picture  was  taken  a  short  time  alter  crack  growth  had  begun. 
As  the  figure  shows,  trying  to  define  the  initial  direction  of  the  crack  is  complicated  by  the  large  deformations,  and 
results  in  ambiguity.  These  complications  are  avoided  by  determining  the  kink  angles  from  the  fracture  surfaces  of 
the  specimens  after  final  fracture. 

As  the  mesh  in  Figure  2  shows,  the  model  is  constructed  without  using  a  “spider  web”  near  the  crack  tip. 
Releasing  nodes  in  front  of  the  crack  tip  allows  for  growth,  but  only  in  a  self  similar  fashion.  A  mixed-mode  crack 
does  not  grow  in  a  self  similar  manner,  however,  but  instead  grows  along  a  curved  trajectory  for  a  short  time  until  it 
orients  itself  with  the  mode  I  direction.  A  precise  analysis  of  this  phenomenon  would  require  a  complicated  analysis 
and  numerous  meshes  since  the  instantaneous  crack  growth  direction  cannot  be  predetermined.  An  alternative 
method  was  used  here,  and  is  shown  in  Figure  4.  For  a  given  specimen,  the  projection  of  the  crack  profile  onto  a 
horizontal  plane  was  measured  from  videotape  images,  and  was  synchronized  with  the  load-displacement  data  so 
that  the  effective  crack  length  and  the  load  were  known  at  different  times  during  the  crack  growth.  Thus  the 
irregularly  shaped  mixed-mode  crack  was  modeled  as  a  flat  mode  I  crack,  which  could  then  be  analyzed  more 
easily.  This  is  justified  by  the  tendency  of  the  crack  to  grow  in  the  mode  I  direction  after  a  small  amount  of  crack 
growth  has  occurred,  and  is  also  suggested  by  the  similar  curve  fits  for  daeffdt  versus  Kj  data  for  all  the  different 
mixities  tested. 

An  added  complication  of  the  crack  growth  analysis  in  rubbery  composites  is  the  sporadic  nature  of  the 
crack  growth.  The  crack  growth  in  this  highly  filled  composite  takes  place  by  nucleation,  growth,  and  subsequent 


6 


coalescence  of  voids  ahead  of  the  crack  tip.  During  this  process,  the  crack  growth  rate  will  be  uneven,  with 
intermediate  periods  of  blunting  and  slow  growth.  Measuring  and  plotting  a^versus  t  in  this  situation  will  give 
growth  rates  that  cannot  be  reconciled  with  any  fracture  parameters  or  remote  loading  conditions.  Therefore 
attempts  to  relate  daejfdt  to  the  remote  loads  through  Kj  will  give  inconsistent  results  unless  the  results  are 
smoothed  before  relating  them  to  crack  tip  parameters  such  as  Kj.  The  investigator  can  use  any  of  several 
techniques  to  smooth  the  crack  size  versus  time  data  before  derivation  of  the  crack  growth  rate  without  influencing 
the  resulting  predictions  significantly.  Techniques  such  as  the  secant  method,  the  incremental  polynomial  method, 
and  the  polynomial  method  have  been  used  in  past  analyses  of  mode  I  cracks.  Using  this  experience  as  a  guide,  we 
have  used  the  polynomial  method.  The  crack  growth  versus  time  data  for  each  trial  is  fitted  with  a  cubic 
polynomial,  which  is  then  differentiated  to  determine  the  crack  velocity,  daefj/dt. 

RESULTS 

Results  for  initiation  of  growth  as  determined  from  the  experimental  results  are  shown  in  Figure  5.  The 
associated  curve  is  an  elliptical  curve  of  the  form: 


K, 


x 2 


\KJCj 


Kn 


=  1  . 


(4) 


Here  the  mode  I  and  mode  II  toughnesses  ( Kjc  and  Kjjc,  respectively)  are  regarded  as  unknown 
parameters  and  are  found  using  a  least-squares  adjustment  of  the  data.147  Equation  (4)  is  regarded  as  linear  with 
respect  to  the  fitting  parameters  ( 1/Kjq )2  and  ( 1/Kjjq )2  and  these  unknowns  are  found  by  minimizing  the  sum  of 
the  square  of  the  residuals.  Assuming  random  variations  in  the  toughness  measurements,  this  procedure  allows  for  a 
determination  of  Kjjc  although  specimens  experiencing  pure  mode  II  conditions  are  never  tested.  As  the  figure 
shows,  the  mode  II  fracture  toughness  is  significantly  lower  than  its  mode  I  counterpart.  This  suggests  that  different 
mechanisms  govern  fracture  as  the  mode-mixity  changes.  However,  there  are  no  sufficient  explanations  of  this 
phenomenon  that  apply  to  a  rubbery  particulate  composite.  Evans  and  Hutchinson  previously  analyzed  cracks  with 
nonplanar  faces,  and  suggested  that  the  crack  face  contact  under  combined  loading  conditions  resulted  in  contact, 
friction,  and  even  frictional  locking.7  *  While  this  explanation  is  feasible  for  many  materials,  this  hypothesis  predicts 
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that  the  pure  mode  II  fracture  toughness  is  higher  than  its  mode  I  counterpart,  which  is  not  the  case  here.  Also, 
during  the  experiment,  the  crack  faces  are  separated  (see  Figure  5),  so  that  crack  face  interaction  does  not  explain 
the  trend  found  here.  Another  near  tip  phenomenon,  plasticity,  has  been  postulated  to  explain  differences  in  fracture 
toughness  with  mode  mixity,  but  again  higher  mode  II  toughness  is  predicted.10  More  research  needs  to  be  done  to 
explain  the  variations  in  fracture  toughness  and  how  this  relates  to  the  micromechanisms  in  the  rubbery  particulate 
composite. 

Although  the  elliptical  failure  locus  is  an  empirical  curve  fit  and  cannot  now  be  linked  to  any  explanations 
based  on  micromechanisms,  the  kink  angles  agreed  approximately  with  several  theoretical  predictions.  All  of  these 
theories  differ  only  in  the  assumed  determinant  of  initial  crack  growth  direction.  Ffyu*  e  6  shows  the  kink  angles 

plotted  against  the  crack  orientation  angle.  The  most  closely  agreeing  theory  is  the  strain  energy  density  theory 

/ 

introduce  in  the  1970’s  by  Sih  and  Chen,5'15  The  quantity  that  determines  the  initial  crack  trajectory  is  the  strain 
energy  density,  which  is  a  function  of  the  two  components  of  the  complex  stress  intensity  factor  K  according  to: 

S  =  anK ]  +  2anKjKn  +  a22K2u  .  (5) 

The  three  coefficients  aj 7,  a] 2,  and  a22  depend  on  9 ;  the  material  properties,  and  the  stress  state.  The 
hypothesis  is  that  the  crack  will  grow  at  an  angle  such  that  the  value  of  S  will  be  maximized.  Letting  cE/dQ  equal 
zero  and  substituting  for  Kj  and  Kjj  from  eqn  (1)  results  in  a  prediction  of  the  kink  angle,  9q ,  for  a  simple  plate 
geometry  that  relates  to  the  crack  orientation  angle  through: 

( k  - 1)  sin(0o  +2f3)-2  sin(2#0  +  2/?)  +  sin(2#0 )  =  0  •  (6) 

Here  k  =  (3-4v)  for  plane  stress  and  (3-v)/(l+v)  for  plane  strain,  where  v  is  Poisson’s  ratio.  Other  theories 
worthy  of  mention  are  also  shown  in  Figure  6;  these  are  the  maximum  tangential  stress  and  maximum  energy 
release  rate  theories.  The  maximum  tangential  stress  theory  assumes  that  the  crack  grows  so  that  oQQ  is  maximized. 
The  maximum  energy  release  rate  theory  analyzes  an  imagined  short  extension  of  a  main  crack  that  grows  from  the 
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main  crack  at  an  angle,10.  The  local  energy  release  rate  for  this  crack  can  be  expressed  in  terms  of  global  factors  (Kj 
and  Kjj  for  the  larger  crack),  and  the  optimal  orientation  of  this  kinked  crack  will  be  that  which  maximizes  the 
strain  energy  release  rate,  G. 

One  conclusion  suggested  by  the  Figure  6  is  that  linear  fracture  mechanics  theories  can  predict,  with  fan- 
accuracy,  the  behavior  of  this  rubbery  particulate  composite.  This  is  true  because  the  overall  material  response  is 
linear  for  the  conditions  considered  here.  This  is  shown  by  the  load  displacement  curves  of  the  cracked  specimens 
(see,  for  example,  Figure  7).  Another  manifestation  of  the  linear  nature  of  the  material  response  is  shown  in 
Figure  8 k16  This  figure  shows  the  measured  and  calculated  crack  opening  displacements  and  also  the  measured  and 
calculated  v  field  fringe  data  for  a  moir£  experiment  performed  on  a  similar  material  (for  simplicity,  mode  I 
conditions  are  shown).  The  analysis  used  linear  fracture  mechanics  concepts  that  are  incorporated  in  the  resulting 
calculated  values  shown  in  the  figure.  The  close  agreement  between  the  experimental  and  the  calculated  fringe 
values  and  crack  opening  displacements  suggests  that  the  linear  fracture  mechanics  assumptions  are  justified. 

Linear  fracture  mechanics  concepts  also  work  well  with  growth  rate  predictions.  Figure  4  shows  a 
growing  mixed-mode  crack  modeled  as  an  equivalent  mode  I  crack.  The  crack  growth  rates,  when  determined  in 
this  way,  can  be  plotted  according  to  a  power  law: 


da 


eff 


dt 


=  ex; 


(7) 


Here  aejf\s  the  effective  crack  length  and  Kj  is  the  mode  I  stress  intensity  factor;  both  are  based  on  the 
model  shown  in  Figure  4b.  Figure  9  shows  the  data  for  all  the  mixed-mode  specimens  plotted  on  a  log-log  graph  of 
daefj/dt  versus  K.  The  first  few  data  points  for  each  specimen  have  been  deleted  so  that  steady  state  conditions  are 
achieved.  The  individual  data  sets  all  have  about  the  same  slope  and  can  be  assessed  as  an  aggregate  so  that  the 
parameters  for  the  crack  growth  law  (when  analyzed  using  a  least-squares  technique)  are  found  to  be  C  =  1.79  x  lO”6 
and  m  -  2.73  for  the  material  studied  here  (assuming  daejy/dt  is  in  mm/s  and  Kj  is  in  kPa  m1/2).  The  pure  mode  I 
data  gives  comparable  coefficients  of  C  =  1.85  x  10*6  and  m  =  2.7 4  We  conclude  that  the  complicated 
geometry  in  Figure  4a  can  be  analyzed  using  the  simpler  geometry  of  Figure  4b  for  die  range  of  mixities  tested. 


V 
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CONCLUSIONS 

Although  rubbery  particulate  composites  have  viscoelastic  properties,  high  elongations,  and  complicated 
failure  mechanisms,  they  can  be  studied,  for  a  given  nominal  strain  rate,  using  the  principles  of  linear  ff  acmre 
mechanics.  When  analyzed  in  this  way,  the  fracture  locus  is  elliptical.  The  initial  crack  growth  angles  match  the 
strain  energy  density  theory  best,  although  for  all  but  the  highest  mode  mixities,  other  theories  gave  nearly  identical 
results.  The  crack  growth  rates  can  be  predicted  using  a  simplified  mode  I  approach.  These  findings  allow  for  a 
simplified  approach  that  circumvents  more  elaborate  and  difficult  analyses  required  by  the  consideration  of 
nonlinear  and  time  dependent  effects.  Future  work  could  include  the  use  of  experimental  methods  to  study  the 
causes  of  fracture  in  these  composites,  especially  the  micromechanisms  near  the  craci  *ip  at  various  mode  mixities. 
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K,  [  kPa  m 1/2] 


kink  angle  [deg] 


^Composites  Part  B:  Engineering 
Timothy  C.  Miller 
Figure  6 

90 
80 
70 
60 
50 
40 
30 
20 
10 

Q_ . . . . . 

0  io  20  30  40  50  60  70  80  9 

crack  orientation  [deg] 

■  max  tensile  stress 
•  max  energy  release  rate 
a  max  strain  energy  density 
x  propellant  experiments 


stress  [kPa] 


‘Composites  Part  B:  Engineering 
Timothy  C.  Miller 
Figure  7 


Composites  Part  B:  Engineering 
Timothy  C.  Miller 
Figure  8 


+  Measured  COD 
-  Calculated  COD 
•  Measured  Fringe  Data 
■  Calculated  Fringe  Data 


